Evidence of quantum criticality in the doped Haldane system Y 2 BaNiOs 
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Experimental bulk susceptibility x(T) and magnetization M(H,T) of the S=l Haldane chain sys- 
tem doped with nonmagnetic impurities, Y2BaNii_ a: Zn I 05 (a; < 0.08), are analyzed. A numerical 
calculation for the low-energy spectrum of non-interacting open segments describes very well exper- 
imental data above 4 K. Below 4 K, we observe power-law behaviors, x(T)~ T~ a and M(H,T) ~ 
T 1 ~ a f a (H/T), with a (< 1) depending on the doping concentration x. This observation suggests 
the appearance of a gapless quantum phase due to a broad distribution of effective couplings between 
the dilution-induced moments. 
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The past two decades have seen a resurgence of in- 
terest devoted to macroscopic quantum phenomena in 
Heisenberg antiferromagnets (HAF). A current issue in 
this field concerns the interplay between quantum spin 
fluctuations and quenched (i.e., time-independent) disor- 
der in one dimension (ID) M. A simple nontrivial model 
to address this issue is the 1D-HAF with spins 5 = 1, 
which has a Haldane gap A ps 0.4J and a short corre- 
lation length £ w 6 in its spin-liquid ground state (J is 
the nearest- neighbor exchange parameter) ^|J^]. In this 
system, the disorder in the form of site depletion leads 
to interesting quantum phenomena due to the creation 
of two effective 5 = 1/2 spins on the opposite edges of a 
segment, near the vacancies ^|||. The inclusion of weak 
magnetic bonds across the spin-vacancies induces bond 
disorder and causes departure from simple finite-size be- 
haviors. In this case, some recent theoretical works, for 
the 5 = 1 chain with variable nonmagnetic doping, found 
a quantum (zero-temperature) phase transition from the 
Haldane phase to a ID random-singlet (RS) phase in 
which effective spins are coupled into singlets over all 
length scales On the experimental side, however, 

this problem has been very little explored from the point 
of view of quantum criticality. While it is known that 
some 5 = 1 diluted compounds can sustain paramag- 
netism down to low temperature, most of the existing 
experimental studies aimed to confirm the existence of 
5=1/2 end-chain states or to reveal the internal struc- 
ture of these end states Q,^] . A noticeable exception is 
the recent observation of a 3D long-range order (LRO) 
caused by nonmagnetic doping in PbNi2V2 0s ||. This 
Haldane system shows, however, substantial interchain 
interactions, so it is almost critical towards the forma- 
tion of a LRO in the absence of intentional disorder. 

In this work, we examine the effects of nonmagnetic 



doping on the bulk magnetic responses in the 5 = 1 
quasi-lD-HAF Y2BaNi05. We chose this well character- 
ized system because of the rather high value of its main 
energy scale, J«280 K, and the very small value of the 
interchain coupling, \ J±/J\ < 10~ 3 fiofl . Furthermore, 
it has a simple orthorhombic crystal structure |Tl[| and 
a small single-ion anisotropy (Z3/J«-0.04, Ej J«-0.01) 
fl0|| . A weak nonmagnetic doping in Y2BaNi05, which is 
realized by substituting Zn 2+ or Mg 2+ ions for the 5=1 
Ni 2+ ions, yields new magnetic states below the Haldane 
gap (AwlOO K) with no sign of LRO down to very low 
temperature |l^-|l4| . Recently, the low-temperature spe- 
cific heat of Y 2 BaNii„ a; Zn a; 05 (x = 0.04) and the ESR 
spectra of a series of Mg-doped compounds have been 
quantitatively explained by an effective model which de- 
scribes the low-energy spectrum of non-interacting open 
5=1 chains 15 16[|. These results and the recent NMR 
experiments in Ref. 8 are strong indications of the exis- 
tence of 5 = 1/2 end-chain states. In the present work, 
we analyze new susceptibility and magnetization data for 
Y 2 BaNii_ :E Zn 2; 05 with 0.04<x<0.08. Above 4 K, the 
simple model of non-interacting segments used in Refs. 
flT5| , |l6[ describes the data very well. However, at very 
low temperature we observe a new regime with divergent 
power law behavior, which is similar to that observed 
in some doped semiconductors fl7|| . Our results support 
the existence of a zero-temperature gapless phase due to 
randomness in the effective bond distribution that is re- 
inforced by the interchain interactions. The low-energy 
effective model for this phase resembles a 2D or 3D RS 
state. 

Figure 1 shows the linear susceptibilities, x(T), mea- 
sured with a commercial SQUID magnetometer for a 
series of polycrystalline samples of Y2BaNii_ x Zn 2 ;05 
(x = 0, 0.04, 0.06 and 0.08) that were prepared through 
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standard solid-state reactions []ll|. The results obtained 
for the nominally pure sample compare well with those in 
the published literature (T3). At low temperatures, x(T) 
is well fitted by the sum of a T-independent term xch a 
Curie law and a thermally activated term with A«100 K 
(see Fig.l). The sum of the activated part and the con- 
stant xo ~ 10~ 4 cm 3 /mol is negligible below 15 K. The 
low-temperature susceptibility (below 10-15 K) for x = 
is mainly due to the Curie contribution that corresponds 
to 1.4 percent of free 5 = 1/2 spins per formula unit. 
This Curie behavior can be interpreted as coming from 
natural crystal defects and a weak excess oxygen. The 
main effect of the intentional nonmagnetic dilution is to 
increase drastically the low-temperature upturn in x{T), 
as is evident from Fig. 1 . This indicates the existence of 
paramagnetic moments whose amount is related to the 
doping concentration. The behavior of x(T) is however 
qualitatively the same, regardless of the level of dilution. 
The data for the doped compounds were first analyzed 
with a modified Curie- Weiss law, x(T) = Xo+C/(T-6), 
between 1.9 and 15 K. A small phenomenological (9«-0.6 
K was required to describe the data. The dependence of 
the Curie constant C as a function of the Zn doping is 
shown in the inset of Fig. 1. At low doping, the observed 
C value is that for two S = 1/2 spins per nonmagnetic 
impurity. At high x, however, the Curie constant should 
be that for one half of a spin-1 per Zn since a short seg- 
ment with an odd (even) number of spins behaves as a 
5 = 1 (5 = 0) 'molecular' entity at low temperature 
and the number of odd chains is half the total number of 
chains. This tendency seems to be observed (see Fig. 1). 
We have also measured the magnetization, M(H,T), for 
applied field up to 5 T. For all samples, M(H,T)-xqH 
shows a non-Brillouin behavior and is smaller than twice 
the Brillouin function for one 5 = 1/2 spin per Zn. 

To go further in our analysis, we have compared the 
experimental data with a theoretical model based on the 
picture of non-interacting segments, as described here- 
after. The energy spectrum of finite segments described 
by the 5=1 Heisenberg Hamiltonian is characterized 
by four low-energy states, a singlet and a triplet ||. It 
corresponds to the existence of a coupling between the 
two 5 = 1/2 edge spins, which may be AF or ferro- 
magnetic (F). The energy difference among these levels 
decreases exponentially with the length of the chain TV, 
as (— l) N e~ N /£, and already for chains of a few lattice 
sites, these four states are separated from the rest of 
the spectrum by an energy of the order of the Haldane 
gap Oj. Therefore at temperature T and magnetic en- 
ergies (J-bH much lower than the Haldane gap, the mag- 
netic and thermal properties of a system composed of 
segments of several lengths are completely described by 
the energy of these four states and the matrix elements of 
the magnetization operator in this reduced subspace. We 
have calculated these quantities using the Density Matrix 
Renormalization Group (DMRG) method. The resulting 



Hamiltonian which describes the low-energy properties 
of a segment including the triplet |1,5 Z ) and the singlet 
state |0) is: 

H eff = E (N) + (Ja(N) + D(3(N)) |0) (0| + D^(N)S^. 
+E 1 (N)(S 2 x ~S 2 y )- f i B Y / H a g av S?, (1) 

where E (N), a(N), (3(N) and j(N) are functions of 
the chain length N (determined from the DMRG data). 
The validity of the last term has been verified explicitly 
by calculating the matrix elements of 5 t + , and 5 t ~~ for all 
chains. The magnetic moment operator is -Vjyi? e //, and 
to compare with experiments in polycrystalline samples, 
we have averaged its expectation value over all possible 
orientations of the crystal. Also, a distribution of chain 
segments corresponding to a random distribution of de- 
fects was assumed [jl6). Figure 2 shows the comparison 
between the measured x(T)T products for T < 12 K 
and the numerical solution. The adjustable parameters 
in the calculations were the doping x and xo- All the 
other parameters were held fixed at the values deduced 
from previous works 0,0): J = 280 K, D/J = -0.039, 
E/J = -0.013, g x = g y = 2.17 and g z = 2.20. The fitted 
values of x (x = 0.041, 0.060 and 0.074) are close to the 
nominal concentrations, and the fitted xo remains in the 
expected range (Xo = 1-0 x 10~ 4 to 1.9 x 10~ 4 cm 3 /Ba- 
mol). Taking into account that uncertainties in the above 
parameters might affect the theoretical curve, the success 
of the anisotropic model of non-interacting segments to 
explain the data above 4 K is remarkable. Note that for 
the cases with nominal concentration of Zn x = 0.06 and 
x = 0.08, a decrease of the order of 5% or less in the ac- 
tual x used, would improve considerably the fitting above 
5 K, but increase the differences between theory and ex- 
periment in the low-temperature part. Instead, because 
of the different rate of change of both curves, it is not pos- 
sible to obtain a good fit for T < 4 K by small changes 
of x or other parameters. A good agreement between 
the experimental magnetization (M (H, T)~xqH) and the 
numerical data above 4 K (not shown) is also achieved 
using the parameters required to describe the suscepti- 
bility data. This noticeable agreement confirms the ex- 
istence of the 5=1/2 end-chain excitations. Note that 
an accurate description of the experimental data needs 
to consider both, a random distribution of defects |l(J 
and the effect of the single-ion anisotropy Jl5|,[l6| , which 
is stronger for shorter segments. 

However, below 4 K, the experimental data deviates 
from the predictions of the above described model. This 
is seen in Fig. 2 for the susceptibility. It is also ob- 
served for our magnetization data M(H, T)~xoH. In 
previous studies of specific heat and electron spin res- 
onance, while excellent agreement was obtained above 
4 K, a discrepancy was also found below this tempera- 
ture and interpreted as the onset of interchain interac- 
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tions not taken into account in this model [|l5|,[L6|. An 
analysis of the susceptibility shows that, in the vicinity 
of 4 K, there is a crossover to a sub-Curie power law 
regime, [\(T) - X a] ~ T~ a with a = 0.83, 0.79 and 
0.76 for x = 0.04, 0.06 and 0.08, respectively (see Fig. 
3). The susceptibility data for the Mg-doped compound 
Y 2 BaNio.95gMgo.o4i05 in Ref. (m) also obey a power- 
law form with a = 0.73 (see Fig. 3). This behav- 
ior, which cannot be reproduced by the model of non- 
interacting segments, is reminiscent of some random ex- 
change Heisenberg systems for which the AF couplings 
between nearest-neighbor S — 1/2 spins are random in 
their magnitude, such as insulating phosphorous-doped 
silicon (Si:P) @. Bhatt and Lee have proposed a the- 
oretical method to explain the properties of Si:P [Esl. 
Within this approach, the spin pairs with the strongest 
AF bonds freeze into inert singlets, leaving behind a new 
ensemble of active spins with a renormalized distribu- 
tion of exchange, and the quantum fluctuations drive 
the whole system into a ground state consisting of ran- 
dom local singlets (RS phase). The behavior of the sus- 
ceptibility, x{T)^T~ a with a < 1, follows from a di- 
vergent power-law distribution of the renormalized ex- 
change, P(J)^J~ a . A scaling of the magnetization has 
been also predicted and experimentally observed [ jl7| , 
M(/f, T)/T 1_Q! = f a {H/T). A similar scaling is, in fact, 
very well obeyed by our data below (but not above) 4 K 
using the a exponents determined from the susceptibil- 
ity. Fig. 3 shows a plot of (M(H,T)-XqH)/T 1 -° 1 versus 
H/T for each of the three samples studied. All the data 
for each sample, for different magnetic fields and tem- 
peratures, lie on a single scaling curve, which is however 
different from the expression of f a (H/T) calculated on 
the basis of the Bhatt and Lee's solution for Si:P |L7|. 
This is illustrated in Fig. 3 for Y2BaNio.g4Zno.06O5. 

These power law behaviors show that the picture of 
the ideal ID anisotropic Heisenberg non-interacting chain 
breaks down below 4 K, and that additional interactions 
play a role. From simple arguments based on perturba- 
tion theory, one expects that the dominant inter-chain 
interaction is a coupling J between S = 1 spins lying 
in nearest-neighbor chains along the b direction, perpen- 
dicular to the chain axis a. A calculation using the cell 
perturbation method gives Jj_b«0.2 K p9| . This inter- 
action should be at least an order of magnitude larger 
than that between two spins on the same atomic chain, 
with a Zn atom in between. As discussed in more detail 
b induces an effective interaction J' be- 
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tween S = 1/2 end-states of neighboring chains along the 
b direction, the magnitude and sign of which depends on 
the detailed position of the defects. The maximum ab- 
solute value of J' is of the order of 10 J±b and it should 
decay exponentially with the distance along the chain be- 
tween the defects. Since the interaction between S — 1/2 
defects lying in different ab planes is expected to be at 
least an order of magnitude smaller, the physics in the 



range 0.1 K< T <4 K seems to correspond to a 2D ar- 
ray of S = 1/2 spins, with random F and AF interactions 
coming from exponentially decaying intra- and interchain 
couplings. 

In agreement with this, the spin glass behavior ob- 
served below 3 K in Ca doped systems Y2_ J ,Ca, / BaNi05 
indicates a dimensionality higher than one in the system 
at low temperature @Q. For a 2D or 3D S = 1/2 
model with random AF interactions, Bhatt and Lee ob- 
tained a power law behavior with a (< 1) depending on 
the doping concentration Jlg| ]. This result is consistent 
with our observations but, as mentioned above, we obtain 
a different scaling function for the magnetization. This 
might be due to the presence of effective F coupling in 
Y2BaNii_ x Zn x 05. Our results can also be compared to 
a recent theoretical work in which the spin-1 chain with 
an AF coupling across nonmagnetic impurities is mapped 
onto a bond disordered spin- 1/2 chain |6). Within this ID 
approach, the Haldane state is found to be stable to weak 
dilution. For sufficiently strong randomness, however, 
the system flows towards a ID RS phase with diverging 
correlation length. In this gapless phase, the form of the 
low energy asymptotic behavior of x(T) is independent 
of the level of the disorder, x(T) ~ l/(Tln 2 T). In the 
Haldane phase, there is a Griffiths gapless region where 
x(T) varies as T~ a with a (< 1) which is dependent 
of the details of the randomness. It is unclear however 
whether a similar physics should be valid or not when 
interchain coupling is turned on pcj . Furthermore, the 
form of the asymptotic behavior of x(T) is not known 
for a 2D or 3D RS phase [Q. If there were no interchain 
coupling J±, the low-energy model for Y2BaNi1_.rZn.EO5 
with 0.04 < x < 0.08 would be the ID RS phase de- 
picted in Ref. |J since the bond across a Zn is expected 
to be much weaker than the average coupling between 
the S = 1/2 edges within a segment |la]. From our 
observations, it is tempting to say that the H/T vari- 
able in the scaling of M(H, T) indicates a T c = critical 
point and that Y2BaNii_ 2; Zn I 05 has a quantum critical 
ground state even for the smallest doping x — 0.04. From 
the theoretical point of view, the detailed explanation of 
the observed power-law behaviors remains open. 

Interestingly, a similar singular behavior has been dis- 
closed for two S = 1 1D-HAF compounds, namely 
AgVP2S6 and NENP, on nominally pure samples [po[ . 
A broad distribution of couplings between the dilution- 
induced moments and quantum fluctuations drives 
Y2BaNii_ :E Zn I 05, and perhaps other real systems as 
well, towards a novel ground state which resembles a ran- 
dom valence bond state. This is in contrast with the LRO 
observed in PbNi2V 2 0s ||. Therefore there is no univer- 
sal behavior under nonmagnetic doping in real 5=1 1D- 
HAF substances. Details of the topology of the magnetic 
interaction lattice are probably crucial in determining the 
way the Haldane phase is destabilized by the quenched 
disorder. 
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FIG. 1. Magnetic susceptibility versus temperature, x(T), 
for Y 2 BaNii_ x Zn ;c 6 (0<x<0.08). The solid line is a fit for 
x=0 with x=Xo+C/(T-6)+A( A/T) 1/2 e- A/T . The sum of the 
first and third terms is shown as a dashed line. Inset: the dop- 
ing dependence of C as described in the text. Data from Ref. 
[12] (open circle) and [14] (square) for Mg-doped Y2BaNiOs 
are included for comparison. Solid and dashed lines show 
the result assuming that one impurity creates two free S=l/2 
spins and one half of a S=l spin, respectively. 
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FIG. 2. xT versus T plot for Y 2 BaNii_ a ,Zn a ,0 5 
(0.04<x<0.08) in the range 1.9-12 K. Errors bars are of the 
order of the sizes of the symbols. The solid lines show the best 
fits with our numerical solution, as described in the text. The 
dashed and dotted lines are for the two free S=l/2 spins pre- 
diction and a numerical calculation neglecting the single-ion 
anisotropy, respectively, for x=0.06. 

FIG. 3. Power-law scaling of the magneti- 
zation, [M(H,T)- X o]/T 1 - a =f a (H/T), for YaBaNi^Zn^Os 
(0.04<x<0.08) in the range 1.9-4 K with data taken every 
0.1K. The values of a are determined from the susceptibili- 
ties. The solid line is the theoretical expression of f a (H/T) 
from Ref. [17] for x=0.06. Inset: Double-logarithmic plot 
of [x(T)-Xo]T versus T for Y 2 BaNii_ a; Zn :l; 5 . Data for 
Y2BaNio.959Mgo.o4i05 in Ref. [14] are shown as open dia- 
monds. 
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